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Abstract 

QSOs allow study of the evolution of the relationship between black holes in galactic nuclei and their host 
galaxies. The black hole mass Mbh can be derived from the widths of the broad emission lines, and the stellar 
velocity dispersion <r* of the host galaxy can be inferred from the narrow emission lines. Results based on [O ill] and 
[O n] line widths indicate that the Mbh — o"* relationship, at redshifts up to z ~ 2, is consistent with no evolution 
or an increase of up to ~ 0.5 dex in A/bh at fixed CO line widths offer an estimate of er» for luminous QSOs at 
high redshifts. The available objects from z ~ 4 to 6 have very massive black holes, Mbh ~ 10 9 ' 5 Mq, but their CO 
line widths suggest much smaller host galaxies than would be expected by the local Mbh — relationship. The 
most massive black holes must continue to reside in comparatively modest galaxies today, because their number 
density inferred from QSO statistics exceeds the present-day abundance of proportionally massive galaxies. 

Key words: galaxies: active galactic nuclei, supermassive black holes 



1. Introduction 

The study of black hole demographics has 
added a new dimension to research involving ac- 
tive galactic nuclei (AGN). This is rooted in two 
developments of recent years. The first is the avail- 
ability of measurements of supermassive black 
holes in nearby galaxies, involving observations of 
stellar and gaseous motions with HST along with 
other techniques (reviews by Kormendy & Geb- 
hardt 2001; Ferrarese & Ford 2004; Combes 2005). 
This has led to the realization that Mbh is closely 
correlated with the luminosity and especially the 
velocity dispersion of the bulge component of the 
host galaxy (Gebhardt et al. 2000; Ferrarese & 
Merritt 2000). 

The local Mbh — a* relationship is given by 



Tremaine et al. (2002) as 

M BH = (lO 8 - 13 M )(<r*/2OO km s" 1 ) 4 ' 02 . (1) 

The rms dispersion of only ~ 0.3 in log Mbh in this 
relationship suggests a fundamental connection be- 
tween the evolution of supermassive black holes 
and their host galaxies. The formation and evolu- 
tion of supermassive black holes has become a focus 
of theoretical study (review by Haiman & Quataert 
2004; Croton 2005, and references therein). 

The second development is the ability to esti- 
mate Mbh in AGN based on evidence that the 
broad emission lines come from gas orbiting the 
black hole at a radius that can be estimated from 
the continuum luminosity (Kaspi et al. 2005, and 
references therein). Central black hole masses in 
AGN can easily be estimated in this fashion, al- 
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lowing demographic studies and providing insight 
into AGN physics. Results below assume a cosmol- 
ogy with H = 70 km s _1 Mpc -1 , Qm = 0.3, and 
n A = 0.7. 



2. Black Holes in QSOs 

Direct measurements of nuclear black holes, 
based on spatially resolved measurements of the 
velocities of stars and gas within the gravitational 
sphere of influence of the hole, are limited to 
nearby objects with negligible look-back times. 
QSOs afford an opportunity to study the Mbh — c* 
relationship as a function of cosmic time. However, 
this requires a measurement of the luminosity or 
velocity dispersion of the host galaxy in addition 
to the black hole mass. The host galaxy luminosity 
and the stellar velocity dispersion are difficult to 
measure directly at high redshift, given the glare 
of the active nucleus. An alternative approach is 
to estimate ct* through the use of some surrogate 
involving emission lines of gas orbiting in the host 
galaxy. 

Nelson & Whittle (1996) found that, on aver- 
age, the width of [O III] tracks the stellar velocity 
dispersion, ct [0 m] = FWHM([0 III])/2.35 « ct*. 
This is supported by the agreement of Mbh and 
C[o in] m nearby AGN with the local Mbh — c* 
relationship (Nelson 2000) . Further support comes 
from the study of Bonning et al. (2005), who found 
overall agreement between <jiq nl ] in low redshift 
QSOs and the value of ct* implied by the measured 
host galaxy luminosity and the Faber- Jackson re- 
lation. The rms scatter is <~ 0.13 in log up in] at 
fixed host luminosity. This suggests that CTp hi] 
may be a useful proxy for ct* in statistical studies, 
although not for individual objects. 

Shields et al. (2003) carried out such a program 
using the narrow [O in] A5007 emission line of 
QSOs. They collected measurements of the widths 
of the broad H/3 line and narrow [O in] line in QSOs 
ranging in redshift from z < 0.1 to z = 3.3. Black 
hole masses were derived from the expression 

M BH = (10 7 - 69 M Q )vl oao L°f, (2) 

where w 30 oo = FWHM(H/3)/3000 km s" 1 and 
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Fig. 1. Deviations from the local Mbh ~ °"* relation- 
ship, A log Mbh, as a function of redshift. Filled squares: 
H/3 - [O III] results from Shields ct al. (2003); small cir- 
cles: H/3 - [O in] and Mg II - [O II] results from Salviander 
et al. (2006); open triangles: [O III] measured from spec- 
tra by Sulentic et al. (2004, see text and Table 1); crosses: 
H/3 - [O III] results from Netzer et al. (2004) and Shemmer 
et al. (2004); large open circles: high redshift CO QSOs 
from Shields et al. (2006), who suggest a recalibration 
downward by ~ 0.5 from the values shown. 

L44 = uL„/(10 44 erg s _1 ), the continuum lumi- 
nosity at 5100 A, based on Kaspi et al. (2000). 
They expressed their results in terms of the de- 
viation of an object's actual Mbh from the value 
M a expected on the basis of Equation 1 using 
a \o in]- We follow their use of A log Mbh = 
log Mbh — log M a . Figure 1 shows the results 
for A log Mbh as a function of redshift. There 
is considerable scatter, but overall the objects 
show little systematic offset from Equation 1 as a 
function of redshift. Shields et al. concluded that 
Mbh at z ~ 2 differs by less than 0.5 dex from the 
present day value expected for a given ct*. How- 
ever, the high redshift objects have larger masses 
(M BH « 10 9 - 6 M Q ) than the low redshift objects 
(~ 10 8 - 3 M Q ), so the evolutionary comparison 
involves some disparity in mass. 

Salviander et al. (2006; see also these proceed- 
ings) used these methods to assess the evolution of 
the Mbh — relationship in the QSOs in Data 
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Release 3 (DR3) of the Sloan Digital Sky Sur- 
vey 1 (SDSS; Abazajian et al. 2005). They used 
H/3 and [O in] for redshifts up to 0.8, and the 
broad Mg n A2800 and narrow [O n] A3727 emis- 
sion lines to reach redshifts up to 1.2. Their results, 
included in Figure 1, nominally show a rise of ~ 
0.5 in A log Mbh from z ~ 0.3 to z ~ 1.0. How- 
ever, they show that two sources of bias contribute 
to this apparent evolution. (1) The scatter in the 
local Mbh — c* relationship causes some galaxies 
to harbor exceptionally large black holes. If these 
are fueled in proportion to their mass, their high 
luminosities cause them to be over-represented in 
the QSO sample. (2) There is a tendency in data of 
marginal signal/noise quality to fail to detect wider 
lines, over-representing narrower [O in] lines in the 
sample. Subtracting these biases, which both act 
to give an increase in A log Mbh with redshift, Sal- 
viander et al. find a residual evolution of 0.3 ± 0.3 
in A log Mbh at redshift unity. 

CO emission lines afford an opportunity to study 
QSO host galaxies at higher redshifts. Shields et 
al. (2006) used ct co = FWHM(CO)/2.35 together 
with Mbh derived from the broad Mg II and 
[C IV] lines to assess the Mbh — c, relationship 
in QSOs up to z — 6.4. They find large black holes 
(~ 10 95 M Q ) that are at least an order of magni- 
tude larger than expected for the CO line width 
(see Figure 1). These giant black holes in the early 
universe evidently reside in comparatively modest 
galaxies, as found for the case of SDSS J1148+5251 
(z = 6.4) by Walter et al. (2005). 

Figure 1 summarizes a variety of measurements 
of A log Mbh as a function of redshift. Included are 
a number of z ~ 2 QSOs from Netzer et al. (2004) 
and Shemmer et al. (2004), as selected by Bon- 
ning et al. (2005). The six radio- loud objects have 
a mean A log Mbh of —0.4. Also shown are nine 
QSOs from the study of Sulentic et al. (2004). For 
these, we have measured the [O in] width from the 
VLT spectra of Sulentic et al., kindly provided by 
P. Marziani (personal communication). We made 
a direct measurement of the FWHM of the A5007 
line for those objects having adequate [O III] in- 
tensity and not having excessive [Fe n] emission, 
keeping nine of 17 objects. We measured the con- 



1 The SDSS Web site is |http://www.sdss.org7| 



tinuum flux at 5100 A rest wavelength from the 
spectra, derived the luminosity (for our cosmol- 
ogy), and calculated Mbh from Equation 2 taking 
FWHM(H/3) from Sulentic et al. The VLT mea- 
surements have a mean A log Mbh = 0.31, or 0.48 
including only the radio-quiet objects. Altogether, 
these results are consistent with the conclusion of 
Shields et al. (2003) that there is relatively little 
evolution of the Mbh — c, relationship between 
z ~ 2 and the present, albeit with considerable 
scatter. 

Tabic 1 

Width of narrow [O III] line in QSOs measured from the 
spectra by Sulentic et al. (2004), along with AfeH derived 
from Equation 2 (see text). Absolute luminosity vL v is 
given at rest wavelength A5100 for our adopted cosmology. 
HE 0003, 0005, 0454, and 2349 are radio-loud (Sulentic ct 
al. 2004). 



QSO z lc 


>g f[0 III] 


log uL v 


log Af B H A 


log M B h 


HE 0003 1.077 


2.41 


46.12 


9.26 


0.68 


HE 0005 1.412 


2.73 


46.25 


9.40 


-0.45 


HE 0048 0.847 


2.41 


45.49 


9.23 


0.64 


HE 0248 1.536 


2.76 


47.49 


9.73 


-0.25 


HE 0331 1.115 


2.34 


46.40 


9.41 


1.11 


HE 0454 0.853 


2.18 


45.83 


8.71 


1.06 


HE 2340 0.922 


2.34 


46.18 


9.07 


0.77 


HE 2349 1.604 


2.83 


46.36 


9.35 


-0.91 


HE 2355 2.382 


2.67 


46.71 


9.77 


0.15 


Average 1.305 


2.52 


46.31 


9.33 


0.31 



3. Does [O in] track er* in AGN? 

We noted above some indications that ctjo in] 
may be a useful surrogate for ct, in a statistical 
sense. Figure 2, based on Bonning et al. (2005), 
compares up in] with ct, in a wide variety of AGN. 
Because of the scatter in CT[o in] at fixed ct, , a wide 
dynamic range in ct, is needed to clarify the overall 
trend. For lower luminosity AGN (Seyfert galax- 
ies), direct measurements of ct, are used. For QSOs, 
in which ct, is difficult to measure, ct* is inferred 
from Mbh and Equation 1. We include here the 
VLT results of Table 1 and, at very low ct, the 
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Fig. 2. Correlation between <T[o hi] an d o", for AGN, from 
Bonning et al. (2005) with the addition of POX 52 (Barth 
et al. 2004) and the VLT data from Table 1. See Bonning 
et al. for details and references. This figure follows Bonning 
et al. in using a slightly different formula for Mbh from 
Equation 2. 

dwarf Seyfert galaxy POX 52 (Barth et al. 2004). 
There is a clear trend of increasing ct[q hi] with u* , 
consistent with the idea that [O in] is a valid, if 
noisy, surrogate. 

Greene & Ho (2005) compare er» with the widths 
of several narrow emission lines in a sample of 
narrow-line Seyfert galaxies from SDSS. They find 
that <7[o ii] agrees in the mean with cr*. However, 
®[o in] exceeds cr* by ~ 0.13 dex unless a correc- 
tion is made for the extended blue wing of the 
[O in] profile. In contrast, Salviander et al. (2006) 
find that cp in] and up n] agree within a few hun- 
dredths of a dex in their sample of SDSS QSOs, 
without any correction for the blue wing. It is im- 
portant to clarify the effect of the blue wing on 
^[o in] m various classes of AGN, and how best to 
correct for it. 



4. Homelessness Amongst the Largest 
Black Holes 

Several authors have noted that the largest black 
hole masses, inferred in the most luminous QSOs 
from broad line widths or the Eddington limit, ex- 
ceed the largest values of Mbh found in the local 
universe, and the implied values of c* by Equa- 
tion 1 exceed the largest c* values in galaxies (Net- 



zer 2003; Wyithe & Loeb 2003; Shields & Gebhardt 
2004). If the black hole masses are correct, this im- 
plies a breakdown of the local Mbh — relation- 
ship at high Mbh- 

McLure & Dunlop (2004), in a study of SDSS 
QSOs, find values of M BH up to ~ 10 10 M Q , but 
they dismiss the largest values as resulting from 
the scatter in deriving A/bh from expressions like 
Equation 2. Thus they conclude that black holes 
in QSOs are consistent with the maximum mass 
Mbh ~ 10 9 5 M© found in giant elliptical galaxies, 
such as M87 (Tremaine et al. 2003). There is indeed 
scatter in the BLR radius as a function of luminos- 
ity (Kaspi et al. 2000), which underlies Equation 1. 
However, we argue here that the large values of 
Mbh in the most luminous quasars are likely real. 
In the QSO sample of Shields et al. (2003), there 
are eight radio-quiet objects with log ^£^(5100) > 
46.6; and of these, five have log Mbh > 9.7, sub- 
stantially exceeding the largest Mbh in local galax- 
ies. Here we arc dealing with a majority of the ob- 
jects of a class selected by a criterion that does not 
involve the scatter in Equation 2. Thus, one can- 
not appeal to scatter to dismiss the large values. 
These six objects have an average log L/L^d ~ 
—0.1, suggesting that Mbh cannot be much less 
than derived from the H/3 width. 

What is the present-day density of these giant 
black holes? From the QSO luminosity function of 
Boyle et al. (2003), we estimate the space density of 
QSOs with log vL u (5W0) > 46.6 to be ~ 6 Gpc" 3 
(comoving) at z = 2. Since they are nearly at the 
Eddington limit, we take their lifetime to be the 
Salpeter e-folding time of 50 million years (effi- 
ciency 0.1c 2 ). Applying the above fraction of 5/8 
and taking an effective QSO epoch of 3 billion years 
(Warren, Hewitt, & Osmer 1994), we find the den- 
sity of relic black holes over 5 billion M Q to be ~ 
10 2 - 3 Gpc" 3 . 

By the local Mbh — cr* and A/bh — Mb u i g e 
relationships (Kormendy & Gebhardt 2001), 
a black hole of 5 billion M Q corresponds to 
cr* wSOOkms- 1 and M b „ig e « 10 12 6 M . The 
largest local giant ellipticals (e.g., M87) have 
cr* w 350 km s _1 (Faber et al. 1997), and the ve- 
locity dispersion function of SDSS galaxies ends 
at ~ 400 km s" 1 (Sheth et al. 2003). Bernardi et 
al. (2006) find only two or three candidate galax- 
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ies in SDSS with er* > 500 km s -1 in a volume 
~ 0.5 Gpc 3 among the objects for which they 
find the evidence for superposition to be weakest. 
The nearest black hole with mass > 10 9 7 M Q 
should be at a distance of ~ 100 Mpc and redshift 
of - 7000 km s" 1 . Wyithe (2006) and Wyithe & 
Loeb (2003) reach a similar conclusion based on the 
assumption that the most luminous QSOs shine 
at the Eddington limit. This distance corresponds 
to the largest cD galaxies, such as NGC 6166 in 
Abell 2199 and NGC 7720 in Abell 2634. Such 
galaxies may be a logical place to look for a 5 bil- 
lion solar mass hole. However, such a black hole 
in NGC 6166 or NGC 7720 would violate Equa- 
tion 1, since these galaxies have central velocity 
dispersions ~ 350 km s _1 (Tonry 1984). 

The space density of galaxy clusters with 
a v > 500 km s _1 (Bahcall et al. 2003) exceeds 
by an order of magnitude our derived density of 
black holes over 5 billion solar masses. Ample 
dark matter halos at this velocity dispersion exist 
in the modern universe, but not individual galax- 
ies. Evidently at this value of 0%, the physics of 
baryon assembly is such that giant black holes can 
form in the early universe but the growth of their 
host galaxies is stunted. Perhaps this involves the 
disruptive effect of the QSO luminosity (Benson 
et al. 2003; Wyithe & Loeb 2003; di Matteo et al. 
2005). The giant holes in the CO quasars discussed 
above appear destined to remain in galaxies of 
comparatively modest proportions. 

We thank K. Gebhardt, J. Greene, and B. 
Wills for useful discussions. EWB acknowledges 
support from a Chateaubriand Fellowship and a 
Pierre and Marie Curie Fellowship. This work was 
supported by Texas Advanced Research Program 
grant 003658-0177-2001; by the National Science 
Foundation under grant AST-0098594; and by 
NASA under grant GO-09498.04-A from the Space 
Telescope Science Institute. 



References 
Abazajian, K., et al. 2005, AJ, 129, 1755 
Bahcall, N. A., et al. 2003, ApJ, 585, 182 



Barth, A. J., Ho, L. C, Rutledge, R. E., & Sargent, 

W. L. W. 2004, ApJ, 607, 90 
Benson, A. J., et al. 2003, ApJ 599, 38 
Bernardi, M. 2006, AJ, in press, astro-ph /0510696 
Bonning, E. W., et al. 2005, ,\|).J. 626. 89 
Boyle, B. J., et al. 2000, MNRAS, 317, 1014 
Combes, F. 2005, |astro-p"h7 0505463 
Croton, D. J. 2005, astro-ph/0512375 
Di Matteo, T., Springel, V., & Hernquist, L. 2005, 

Nature, 433, 604 
Faber, S. M., et al. 1997, AJ, 114, 1771 
Ferrarese, L. & Ford, H.. 2004, Space Science Rev. 

116, 523 

Ferrarese, L. & Merritt, D. 2000, ApJL, 539, L9 
Gebhardt, K., et al. 2000, ApJL, 539, L13 
Greene, J. E., & Ho, L. C. 2005, ApJ 627,721 
Haiman, Z., & Quataert, E. 2004, in Supermassive 

Black Holes in the Distant Universe, ed. A. J. 

Barger, Astrophysics and Space Science Library 

(Dordrecht: Kluwer), Vol. 308, p. 147 
Kaspi, S., et al. 2000, ApJ, 533, 631 
Kaspi, S., et al. 2005. ApJ, 629, 61 
Kormendy, J. & Gebhardt, K. 2001, AIP 

Conf. Proc. 586: 20th Texas Symposium on Rel- 

ativistic Astrophysics, 586, 363 
McLure, R. J., & Dunlop, J. S. 2004, MNRAS, 352, 

1390 

Nelson, C. H., & Whittle, M. 1996, ApJS, 99, 67 
Nelson, C. H. 2000, ApJ 544, L91 
Netzer, H. 2003, ApJL, 583, L5 
Netzer, H., et al. 2004, ApJ, 614, 558 
Salviander, S., Shields, G. A., Gebhardt, K., & 

Bonning, E. W. 2006, submitted. 
Shemmer, O., et al. 2004, ApJ, 614, 547 
Sheth, R. K. 2003, ApJ, 594, 225 
Shields, G. A., & Gebhardt, K., 2004, Bull. A AS, 

204.6002S 

Shields, G. A., et al. 2003, ApJ, 583, 124 
Shields, G. A., Menezes, K. L., Massart, C. 

A., & Vanden Bout, P. 2006, ApJ, in press, 

astro-ph/0512418 
Sulentic, J. W., et al. 2004. AstrAp, 423, 121 
Tonry, J. L. 1984, ApJ, 279, 13 
Tremaine, S., et al. 2002, ApJ, 574, 740 
Walter, F., et al. 2004, ApJL, 615, L17 
Warren, S. J., Hewitt, P. C, & Osmer, P. S. 1994, 

ApJ, 421, 412 
Wyithe, J. S. B. 2006, MNRAS, 365, 1082 



5 



Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 595, 614 



G 



